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Complete solutions are given in a closed analytic form for unambiguous discrimination of three 
general pure states with real mutual inner products. For this purpose, we first establish some general 
results on unambiguous discrimination of n linearly independent pure states. The uniqueness of 
solution is proved. The condition under which the problem is reduced to an {n — l)-state problem is 
clarified. After giving the solution for three pure states with real mutual inner products, we examine 
some difficulties in extending our method to the case of complex inner products. There is a class of 
set of three pure states with complex inner products for which we obtain an analytical solution. 



o 

Oh' 
(D 

o 
m 

Oh' 



> 



in 

o 
o 



X 



I. INTRODUCTION 

In quantum mechanics, there is no way to distinguish 
nonorthogonal states perfectly. This is because quantum 
measurement is statistical in nature, and it generally de- 
stroys the state of the system to be measured. Consider- 
able works have been done to find optimal ways to distin- 
guish quantum states in various situations [l|-y|- These 
studies are stimulated by needs of quantum communica- 
tion and quantum cryptography where the discrimination 
of quantum states is one of the key issues. 

In the quantum state discrimination problem, two set- 
tings have been commonly investigated. In minimum- 
error discrimination [l| , the discrimination success prob- 
ability is maximized without any constraint on the prob- 
ability of erroneous results. In unambiguous discrimina- 
tion, however, the success probability is maximized under 
the condition that measurement should not produce erro- 
neous results. This is possible by allowing an inconclusive 
result: "I don't know." Ivanovic [j|, Dieks Q, and Peres 
[6| found the optimal solution for unambiguous discrimi- 
nation of two pure states with equal occurrence probabil- 
ities. Later the case of general occurrence probabilities 
was solved in Ref. [7[. Unambiguous discrimination was 
experimentally demonstrated in Ref. [8[. A scheme that 
interpolates the two discrimination settings has also been 
proposed, and its solution for two general pure states was 
given in Rcfs. ^^. 

For unambiguous discrimination among more than two 
pure states, it is not easy to obtain analytical solutions, 
though a great deal of works on general theories have 
been reported: For example see Refs. pT] - [l9| . The same 
thing is true for unambiguous discrimination for mixed 
states [23-i3- 

The case of three pure states has also been consid- 
ered as an application by some articles dealing with gen- 
eral theories. Even for three pure states, to the best 
of our knowledge, complete analytical results are known 
only in limited special cases: (i) The case of symmetric 
states with equal occurrence probabilities [12[, in which 
the states are generated by a single unitary operator; 
(ii) the case in which two of the three real mutual inner 
products are equal and equal occurrence probabilities are 
assumed [ij ; and (hi) the case of three pure states with 



general occurrence probabilities and one of the mutual 
inner products being |19j . 

The purpose of this article is to give complete ana- 
lytical solutions for unambiguous discrimination of three 
pure states with real mutual inner products and general 
occurrence probabilities. In unambiguous discrimination, 
the optimal strategy may produce a vanishing probabil- 
ity of identifying the input state with some of the states. 
This issue will be thoroughly analyzed. 

In Sec. [ni we first establish some general results on 
unambiguous discrinhnation of n pure states. We formu- 
late the problem as semidefinite programming. It turns 
out very useful to consider a problem in which some con- 
straints in the problem are relaxed. This relaxed problem 
is used to clarify the conditions under which some identi- 
fication probabilities vanish, and the problem is reduced 
to a certain {n — l)-state problem. Then, complete an- 
alytical solutions of three pure states with real mutual 
inner products are presented in Sec. IIIII In Sec. IIV| we 
discuss some difficulties in extending our method to the 
case of general complex inner products. In the appendix, 
we show that solutions of unambiguous discrimination for 
linearly independent n pure states are generally unique. 



II. GENERAL CONSIDERATION 
A. Formulation of problem 

Suppose that a state is drawn from a set of known n 
pure states \4>i) with probabilities rji. We assume that 
the n states are linearly independent and the occurrence 
probabilities rji are nonzero. The task is to unambigu- 
ously identify the given state with one in the set of n 
states: An erroneous identification is not allowed, but an 
inconclusive result ( "I don't know" ) is possible. 

Measurement is described by a positive operator- 
valued measure (POVM), and it consists of n -I- 1 ele- 
ments {Ei,E2t . ■ , En, i??}, each element being a posi- 
tive semidefinite operator on the n-dimensional space V 
spanned by \4>i). When the outcome is i {i = 1, . . . , n), 
the given input state is identified with the state | (/>i ) ■ 
The element E-t corresponds to the inconclusive result. 



The POVM satisfies the following completeness relation: 

71 

e. + J2e, = 1v. (1) 



since | (j)i ) are linearly independent. 

Thus, the unambiguous discrimination problem of the 
ensemble {rji, \ (j)i )}"^i is formulated as semidefinite pro- 
gramming: The problem is to maximize 



The no-error conditions are written as 

(0,|i?.|0,)=O (*^j), (2) 

which require that the support of the positive semidef- 
inite operator Ei be the one-dimensional subspace that 
is orthogonally complement to the (n — l)-dimensional 
subspace spanned by the states {\^±)}j{^i)- Therefore, 
Ei should take the following form 



Ei = Xi\(j>i){(t)i I, 



(3) 



where Xi is a non- negative constant to be determined, 
and the state | i^i ) is orthogonal to any state | (f>j ) for 
j ^ i. For convenience we normalize the states | (j)i ) such 
that 

{4'^\4'J)=S^,. (4) 

Since the states | (pi ) are linearly independent, those 
states I (pi ) are uniquely given by 






(5) 



where 



N,: 



■J'l \V] h 



is the Gram matrix of the states to be discriminated. 

The coefficient Xi represents the conditional probabil- 
ity of successfully identifying an input state given that 
the input state is | 0^). In terms of Xi, the discrimina- 
tion success probability is given by 



p = X! ''''^ "^^ i-^''i "^^o = X! ^^ 



(6) 



Now let us examine what condition the positivity of 
E-r further imposes on the coefRcients cc,;. By the com- 
pleteness relation ([T|), the positivity of E-; is expressed 
as 



Iv -Y,x,\^,){4),\ >0. 



(7) 



This condition is equivalent to the positivity of the nx n 
matrix M defined by 

= Nij - XiSij, (8) 



i=l 



subject to conditions given by 



Xi>0 (i = 1,2, . . . ,n), 
N-X >0, 



(9a) 



(9b) 
(9c) 



where X = diag(a;i, X2, ■ ■ ■ ,Xn) and variables are Xi. We 
can show that solutions Xi are uniquely determined; the 
optimal measurement is unique. The proof is presented 
in the appendix. 



B. Relaxed problem and reduction of the number 
of states 



In unambiguous discrimination, the optimal strategy 
may produce a vanishing probability of identifying the 
input state with some of the states \(pi)'. Some of Xi may 
be zero in solutions of the problem (|9]). To see when this 
happens, it is convenient to consider the problem where 
the conditions (|9b|) . Xi > 0, are omitted. To distinguish 
the two problems, this relaxed form of the problem is 
called a "relaxed n-state problem" from now on, and the 
original one (jH]) is simply called an "n-state problem." 

In the relaxed n-state problem, the task is to maximize 



X!^'^^ 



subject to conditions given by 

N-X^ > 0, 



(10a) 



(10b) 



where X^ = diag(a;5^, Xj , . . . , x^). We use a superscript 
R to represent quantities in the relaxed problem. 

There are important relations between solutions of the 
two problems. Let a;^ = {x^, x^, . . . , x^) be a solution 
of the relaxed n-state problem. Then we can show the 
following properties: 

(I) If xf > for all i, then x^ is the solution of the 

n-state problem. 

(II) If xf- < for i in some nonempty index set S, then 

the solution of the n-state problem has vanishing 
components, Xi =0, for at least one i in S. 

Property (I) is evident. To show property (II), we con- 
sider a linear interpolation between the two solutions as 



x{k) = (1 — k)x + KX^ , < K < 1, 



(11) 




n-state problem 



relaxed 
n-state problem 



FIG. 1: (Color online) Linear interpolation of the solutions of 
the n-state problem and the relaxed n-state problem. Each 
component of x{k) is plotted vs. k. There exists < kq < 1 
such that Xk{Ko) > for all k and Xi{iio) — for a certain i G 
S. By the uniqueness of the solution of the n-state problem, 
Ko is actually 0. 



where x is the solution of the n-statc problcni and x^ is 
a solution of the relaxed n-state problem. We find that 
x{k) satisfies the condition (Pc]) since 

N - X{n) = (1 - n){N - X) + k{N - X^') > 0. (12) 

The success probability corresponding to x(k,) is written 
as 



P{k) = (1 - K)pmax + f^Pmax; 



(13) 



where Pmax and p^Ji^x a-re the maximum success probabil- 
ities of the n-state and relaxed n-state problems, respec- 
tively. Note that p^ax ^ Pmax by definition. From this 
we note that p{k) is no less than Pmax for < k < 1. This 
implies that x{k) is a solution of the n-state problem if 
its all components Xi^n) are nonnegative. 

Obviously, ii(/v) is continuous with respect to n. We 
also know that Xi > for all i and that xf < ioT i £ S, 
whereas ccf > otherwise. Using the intermediate value 
theorem, we observe that there exists < kq < 1 such 
that a;fc(Ko) > for all k and Xi(Ko) ~ for a certain 
i & S (see Fig. [1]). Therefore, x{ko) is a solution for 
the n-state problem. Remember that the solution of the 
n-state problem is unique, implying kq ~ 0. Thus, we 
can conclude that the solution of the n-state problem 
has vanishing components, Xi — 0, for at least one i in 
S. This completes the proof of property (II). 

Next, we show that if some of x^ are negative in a solu- 
tion of the relaxed n-state problem, the n-state problem 
can be reduced to a certain (n— l)-state problem. Let us 
first assume that only one component is negative. With- 
out loss of generality, we assume that x^ < and the oth- 
ers are non-negative. Property (II) then states that the 
component Xn is in the solution of the n-state problem. 



This solution can be obtained by requiring an additional 
constraint, x„ = 0, to the original n-state problem. With 
this additional constraint, the condition ^ becomes 



ri-l 

Let us define the following two sets of n — 1 states: 



^/{^^\Q\W^ 



(z = l,2,.. 



1), 



(14) 



(15) 



I </):;) = 7(07101^1 </),: > (*-l,2,...,n-l), (16) 
where 

g = iy-|0„)((/.„|, (17) 

is the projector onto the (n— l)-dimensional subspacc V' 
spanned by the n — 1 states {| (pi )}"^i- 

It is readily verified that the two sets of states are 
reciprocal in the sense that ( (pi \ 4''j ) ~ ^ij ■ Furthermore, 
the condition p4)) can be expressed as 



lv'-J2x',\^',){4>',\>0, 



i=l 



where x' is defined as 



\Q\' 



(18) 



(19) 



We will shortly see that x'^ are variables in the reduced 
(n — l)-state problem. The success probability of the n- 
state problem with the additional constraint, a;„ = 0, is 
rewritten as 



n-l 



where 77- is defined as 




/n-\ 



51 ^i^;,- 



T\i{(i>i \Q\^i 



En— 1 / 



\Q\' 



(20) 



(21) 



Note that if^ > and X^i^i v'i ~ 1- The ri[ turn out to be 
the occurrence probabilities in the (n— l)-state problem. 
From Eq. (PU)) and the condition p^ . we observe that 
the n-state problem with the additional constraint Xn = 
is equivalent to the (n— l)-state problem for the ensemble 
{r][^ l</'i)}i!=i • The maximum success probabilities and 
the variables Xi of the two problems are related in the 
following way: 

p'£lA{^^^4'^)ru] 

= (E^^('^^IQI^^) Pmax'^ [HAmiZl] , (22) 



\k=l 



Xi = {(f>i |Q|0i)a;-. 



(23) 



It should be emphasized that the extra condition x„ = 
does not imply that the state | </>„ ) can completely be 
neglected. This is because the no-error conditions ([2]) 
should be respected for all the states including | 0„ ) . The 
states to be discriminated in the reduced {n — l)-state 
problem are not {| (fn )Yl^i but {| <!)[ )}"^i^, which are the 
normalized states with their | 0„ ) component subtracted. 
As can be seen from (PT|) . the occurrence probabilities are 
also modified by the factor of {(t>i \Q\4)i). 

So far, we assumed that only one component of x is 
negative. If more than one component is negative, we 
do not know which Xi {i £ S) vanishes in the n-state 
problem. Therefore, for each i in S, we need to solve the 
corresponding {n — l)-state problem with Xi being set 0. 
From Eq. (|22p . we obtain a finite number of candidates 
for the maximum success probability of the n-state prob- 
lem. Evidently, the one with the largest probability is 
the solution of the n-state problem. 

We can show that the number of negative components 
of x^ is n — 1 at most. This can be seen in the following 
way: The Gram matrix N is strictly positive definite 
since the states {| (pi )}"^i are linearly independent. This 
implies that the condition N — X^ > is satisfied for 
a set of sufficiently small but positive xf^, which leads 
to a positive success probability. Therefore, a;^ with all 
components negative cannot be a solution of the relaxed 
n-state problem. 



C. Dual problem of the relaxed problem 

In our method, we first solve the relaxed n-state prob- 
lem defined in Eq. ([TU)) . If all zf are non- negative, they 
give the solution of the n-state problem. As shown in the 
preceding section, if some of xf are negative, the prob- 
lem is reduced to a certain (n — l)-state problem. How, 
then, do we solve the relaxed n-state problem? We find 
it convenient to consider the dual problem of the relaxed 
n-state problem. 

The task in the dual problem is to minimize 



d = tr YN, 
subject to the conditions given by 
>">0, Yu^^]^, 



(24a) 



(24b) 



where the variable is an n x n Hermitian matrix Y . It 
is readily shown that the minimum dmin gives an up- 
per bound for the maximum pj^j^x ^^ the relaxed n-state 
problem (fTO)) . 



i—1 n— 1 

= tr YX^ = ti Y[N - (N - X^)] 
<trr7V = d^in, (25) 



where the inequality in the last line follows from the in- 
equality try(iV - X^) > 0. Therefore, the set of condi- 
tions given by 



{N - X^)Y = 0, 



(26) 



together with the conditions (|10b[) and (|24b[) is a suffi- 
cient condition for (i,„i,i to be equal to pj^ax- I^i the next 
section, for the n = 2 and n — 3 cases, we first solve the 
dual problem (|24p and then we explicitly verify that this 
set of conditions are satisfied. According to the general 
theory of semidefinite programming |28| , it is guaranteed 
that the solution of the n-state problem can always be 
obtained in this way. 

The variable matrix Y in the dual problem is con- 
strained by the condition ()24bl) . Any positive semidef- 
inite matrix can be expressed as the Gram matrix of 
some set of vectors. Therefore, we can parametrize Y 
satisfying the condition (j24b[) as 



Y 

-L 1.'] 



> Vi I V] , 



(27) 



where | yi ) are arbitrary vectors with | y; ) — ^/rjl- This 
convenient parametrization will be used in the following 
sections. 



THREE PURE STATES WITH REAL 
INNER PRODUCTS 



III. 



In our scheme, we solve the dual problem ([24| to ob- 
tain the solution of the relaxed problem ([TU]). If all xf 
are non-negative, they give the solution of the original 
nonrelaxed problem (|9]). Otherwise, the problem is re- 
duced to the one for n, — 1 states. To demonstrate how 
this scheme works, we first solve the problem of two pure 
states by our scheme. This provides a helpful guideline 
for the three-state problem which will be discussed next. 



A. Two pure states 

The task of the dual problem is to minimize 



subject to 



d = tr YN, 



Y>Q, Yu = rj. 



(28) 



(29) 



where variable y is a 2 x 2 Hermitian matrix. By chang- 
ing the phase of the state | 0i ) , we can assume that the 
mutual inner product is real and negative, ( 01 | 02 ) = 
— |( 01 I 02 ) I ■ Using the parametrization ([27| , we can 
write d as 



rf= 1 - ((yi |2;2 ) + (2/2 1 yi)) 1(01 102)1 



(30) 



Applying the triangle inequality and Schwarz's inequality 
to the second term of Eq. ([5(I)) . we have 

|(2/i|2/2) + (2/2|yi)| <2|(yi|2/2)| 



< 2v/(yi 1 2/1 ) v/( y2 1 2/2 ) 
= 2y/mv2, 

where {yi\yi) = 77^ is used. By this inequality, we obtain 
the minimum of d to be 



rfmin = 1 - 2^771772 I ( (/'I I 02 ) I 



(31) 



and the vectors | yi ) attaining the minimum can be taken 
to be the following one-dimensional vectors: 



I yi ) = Vm, I y2 

Thus, Y is given by 



Y = , . 

with /3 defined to be 

/3 = 



7/2 ' 



(32) 



The attainability condition is given by Eq. (|26p . Sub- 
stituting Y in Eq.p6|) and multiplying this equation by 
/3 on the right, we obtain 



[N -X^)I3 = 0, 



(33) 



since 0^ (3 = 1. By using the matrix representation, 
Eq. ((33|) is rewritten as 



1 — X 



1 

l\(t>2\ 



l(0l|02)| 



From this equation, we immediately find that x^ and x^ 
are given by 



2-1 



X2 



l(0l|02)| 
l(0l|02)| 



The xf must further satisfy Eq. (|10bp : The eigenvalues 
of A^ — X^ should be non-negative. The characteristic 
equation of A — X^ is given by 



A A- 



Vvim 



0, 



which clearly shows that the two eigenvalues are non- 
negative. Therefore, xf is the solution of the relaxed 
two-state problem. 



Thus, if x^ are all non-negative, they give the optimal 
solution of the two-state problem Xi as 



(34) 
(35) 



and the maximum success probability Pmax is given by 



Xi = 


= 1 - 


{fM.l<^.)\. 


X2 = 


= 1- 


^JMMl. 



Pmax = 1 - 2^7/1772 \{ (j>l \ (t>2)\ ■ 



(36) 



If one of xf is negative, however, omitting one of the 
states to be discriminated is optimal. Let us first assume 
x^ < 0, which happens when the occurrence probabilities 
and the inner product satisfy the following condition: 



<l('^l|02)|. 



(37) 



We then need to search for the optimal solution with 
X2 = 0. According to Eqs. (fTS)) and ((2T|) . the two-state 
problem with the additional constraint X2 = is equiva- 
lent to the one-state problem for the following state | (f>[ ) 
with occurrence probability ?][■. 



Oh 



vWW^i 



Vi 



vi{<Pi \Q\(f>i 
??i(0i \Q\(f>i 



1, (38) 



where Q is the projector defined by Q = 1 — | 02 ) ( '/'2 I ■ 
The maximum success probabilities and the variables Xi 
and x'^ of the two problems are related as follows: 

p'^Uim, I 0. )}Li] = ^i(0i 101 01 )pilL[{ri[, I 0'i )}], 
xi = {(f>i \Q\ 01 )x[, 

where x[ is the solution of the one-state problem. 

The one-state problem is trivial: PmlxiWn I 0i )}] = 1 
and x'l = 1, since there is only one state to be discrim- 
inated. Thus, the maximum success probability Pmax is 
obtained to be 



Pn 



r/i 1-1(01 102 )f 



and Xi {i = I, 2) is given by 

x, = l-\{ 
a;2 = 0. 



(39) 



(40a) 
(40b) 



We assumed that x^ < 0. For the case of x^ < 0, 
it is clear that this problem becomes a similar two-state 
problem with the additional constraint xi ~ 0. 

The maximum success probability and the optimal so- 
lution Xi arc summarized in the following way: 



(Z) J^< 1(01 102)1 
771 

Pmax = 77l(l-|(0l|02)|'), (41) 



X2 = 0. 



(42a) 
(42b) 



(**) |(</'l|'/'2>|<,/^< 



1 



\<?^l 192^1 



(43) 



(iii) 



a^l = l-^;^l(0l|'/'2)|, 


(44a) 


^2 = 1-^1 1(01 |</'2)|. 


(44b) 


1 . h2 

l(0l|02)| V'?l 




Pmax = ??2 (l- 1(01 |02)|^j , 


(45) 


Xi =0, 


(46a) 


X2 = l-|(0l|02)|'. 


(46b) 



This reproduces the well-known result of the unam- 
biguous discrimination problem between two pure states. 
These results will be used in the subsequent section. 



B. Three pure states 

We consider the unambiguous discrimination problem 
between three pure states, |0i), |02), and |03), with 
occurrence probabilities r/i, 772, and 773, respectively. 

The inner products between the states depend on 
phases of the states. However, F defined by 



! 01 I 02 )( 02 I 03 )( 03 I 01 



(47) 



is invariant under any phase change of each state. A little 
thought shows that the phases of the three mutual inner 
products can be chosen to be the same by appropriate 
phase changes of the states. 

Here, we assume that F is real, implying that the three 
mutual inner products can be chosen to be real. We 
separately consider two cases of F < and F > 0. The 
case of complex F will be discussed in the next section. 



1. Case of r<0 

In the case of F < 0, the phases of the mutual inner 
products can be chosen to be real and negative. As we 



will see, we can proceed along the same lines as in the 
two-state problem. The d is given by 

d=tTYN 

=i-((2/i|y2) + (2/2|yi))|(0i|02)| 

- ((2/2|y3) +( 2/3 I y2)) 1(02 I 03)1 

-((2/3|yi) + (2/i|y3))|(03|0i)|. (48) 

Using the triangle inequality and Schwarz's inequality as 
in the two-state case, we immediately obtain the mini- 
mum of d to be 



1-27^^1(01 I 02)1 -27^^1(02 I 03)1 



2V^^I(03l0i)l 



(49) 



and the vectors | yt ) that attain this lower bound are the 
following one-dimensional vectors: 

I yi > = Vm, I y2 ) = Vm, I ys ) = Vm- 

Thus, Y can be expressed as 

Y = f3f3\ (50) 

where /3 is defined by 

/3 = 

The attainability condition (p6|) again takes the form, 
{N - X^) f3 = 0. (51) 

since /3^/3 = 1. In the matrix representation, Eq. ([STjl is 
written as 




1-xf -1(01 I 02)1 -1(03 I 01 )| 

|(0l|02)| 1-xf -1(02103)1 

|(03|0l)| -1(02103)1 1-xf' 




==0. 



Solving this equation, we obtain xf , x^ , and X3 to be 



x- = 1 - 



(V^|(0l|02)|+^/^|(03|0l)|), 



(V^l(0l|02)|+V^l(02|03>l), 

(^A7^l(03|0l)| + ^A^|(02|03)|). 



Now we check if x^ satisfy the condition ()10b|) . that 
is, the positivity of A^ — X^. The attainability condition 
([5T|) shows that one of the eigenvalues of A^ — X^ is 0. 
The remaining two eigenvalues are determined by the 
following quadratic equation: 

A^ - oA -^ 6 = 0, 



where the coefficients a and b are given by 

a = I ( 01 I (/>2 ) I + I ( 02 I (^3 ) I 



V3 + Vl 



Vnva 



y/mm 



\{h\^i)\- 



^_ |(02|03)(03|0l)| ^ |(0l|02)(03|0l)| 
V'71'72 V'?2f?3 

, |(0l|02)(02|03)| _ 

The two roots, the eigenvalues of the real symmetric ma- 
trix, must be real, and the quadratic formula shows that 
they are both non-negative since a > and 6 > 0. Thus, 
the positivity of TV — X^ is satisfied, and xf' is therefore 
the solution of the relaxed problem. 

If xf' are all non-negative, xf is also the optimal solu- 
tion of the original three-state problem. 



xi = 1 — 



Xa = 1 



X3 ^ 1 



1 



(V^|(0l|02)| + V^|(03|01>|), 



(52a) 



(V^l(0l|02>|+V%l(02|03>|), 

(52b) 

m\{<l>3\4>i}\ + Vm\{<t>2\h)\), (52c) 



and the optimal success probability Pmax is given by dn 
as follows: 



Pn 



1 - 2V^^|( 01 I 02 )| - 2V^^|( 02 I 03 )| 

-2V^7^|(03|0i)|. (53) 



If some of xf are negative, the problem can be re- 
duced to a two-state problem. The two-state ensemble 
WiJ I 4>i )}?=i to be considered is determined by general 
formulas of Eqs. ([T5]) and (PT|) . The maximum success 
probability and the solution Xi of the original problem 
can be obtained by Eqs. ([^^ and ([231) in terms of the 
well-known solutions of two-state problems. 



2. Case o/ r > 

In the case of F > 0, the phases of the mutual inner 
products can be chosen to be real and positive. We in- 
troduce three positive real numbers ai, 02, and as as 



"1 



^2 



"3 



^(01 


1 02 ) ( 03 


01 ) 


( 02 1 03 > 


^01 


1 02 > ( 02 


03) 


(03|0l) 


/(02 


1 03 ) ( 03 


01 ) 



; 01 I 02 ) 



(54) 
(55) 
(56) 



By using ai, the inner products can be expressed as fol- 
lows: 



( 01 I 02 ; = aia2, ( 02 I 03 / = 0:20:3, 

The d can then be rewritten as 
3 

3 

= l-F^a,aj(yj|2/j) 



asQfi. 



1 - Y^ r],al 

t 



3 

E 

1=1 



OiWt 



(57) 



Now, we must minimize d. This means that we must 

X 3 

minimize 



which is the norm of the sum of 

three vectors ai | j/i ) , a2 1 2/2 ) , and as 1 2/3 ) . Here | yi ) is 
any vector with the norm y^. The minimum value of 



Ei=ia«l2^i 



is determined depending on whether the 



three lengths Oi^fql satisfy the following set of relations: 



(58a) 
(58b) 



ai^/rjl < a2v^ + "3v^> 
0:2 V^ < oi^/ql + o.z^Jrfy, 
Ozy/rji < ai^/rf[ + 02^/^. 



Be) 



These relations are the condition for the three vectors 
Oi\yi) to form a triangle, and we call it the triangle con- 
dition hereafter. 

First, we assume that the triangle condition 



satisfied. 



Ei=iai|yi 



Ml) is 
is 0, 



The minimum value of 

achieved by one-dimensional complex vector | j/j ) = yi 
with I?;; I = ^/rji such that aij/i -I- a22/2 + o^y^ = 0. The 
minimum value of d is then given by 



^min -L / ^ 



ViO, 



(59) 



i=l 



Writing 




we have Y = y*y^ . Thus, Eq. (pS)) is rewritten as 

{N - X«) y* = 0, 

from which we obtain the following three equations for 
x^ , x^, and x^: 



„R\ „,* 



1 - Xi) 2/1 + ai"22/2 + oia3y^ = 0, 



R'\ „,* 



a2aiyl + (1 - x^) 2/2 + "2032/3 = 



RA „,* 



a^aiyl + a3a22/2 + (1 " X3) 2/3=0 



(60) 
(61) 
(62) 



Now, remember that y, are chosen in such a way that 
aiiji + 022/2 + a3?/3 = 0. Thus, the first equation is given 

by 



1 



2\ * 



0. 



Since yl ^ 0, wc obtain x^ — l — af. Similarly, we obtain 



x^ and x^. 



xf = 1 -al, 
xf = l-al 



We can easily see the positivity oi N — X^, since N - 
X^ is expressed as ctot^ , where a is defined by 



a = 




Thus, xf we obtained is the optimal solution of the re- 
laxed three-state problem. 

If xf- arc all non- negative, they are the optimal solution 
Xi of the original three state problem. 



a:i = 1 — ai, 
X2 = I - al, 
xa = 1 — Qq. 



(63a) 
(63b) 
(63c) 



and the success probability Pmax is given by dann as fol- 
lows: 



Pn 



= 1 - {r]ial + ri2al + 'q^al) 



(64) 



If some of xf turn out to be negative, we can solve the 
problem by reducing it to the two-state problem as de- 
scribed before. 

Next, we consider the case in which the triangle 
relation (|58p is not satisfied. Comparing the norms 
of tti I j/i ) , a2 1 2/2 ) , and as 1 2/3 ) , we realize that the 



norm 



E,=iai|y*) 



is minimized by subtracting the 
two smaller norms from the largest one. Suppose that 
o^iy^ — oi2\/rJ2i oi^^/rj^. Then, the minimum value dmin 
is given by 

3 

4=1 

= 1 - 2 ^771 772 a 10:2 + 2^772 ry3a2a3 - 2^773 ?yia3ai, 

(65) 

and the vectors | yi ) that attain this minimum value are 
the following one-dimensional vectors: 

\yi) = V^: \y2) = -Vm, \y3) = -Vv3- 

Thus, Y is given by /3/3' with f3 defined to be 



/3 




The rest of the argument proceeds in a quite similar 
way to that in the case of F < 0. Since Y = (3(3\ the 
attainability condition Eq. ([26]) can be rewritten as 



{N - X^) f3^0, 
with its matrix representation given by 



(66) 



1 



xf 



aia2 oiias 



a2ai 
asai 



1 



x 



R 



^203 



03012 1 




= 0. 



Solving this equation, we obtain xf as follows: 



xf^l 



X2 



Ol 



02 



(a2V^ + "3 073) 



=: 1- 



03 



(O] 
(O] 



m 



a2\/rJ2) ■ 



The positivity of A^ — X^ can be verified in the follow- 
ing way: One of the eigenvalues of iV — X^ is 0, which 
can be seen from Eq. (j66p . The remaining two eigen- 
values are given by the roots of the following quadratic 
equation: 



A^ - aA + 6 = 0, 



where 



Ol 



(a2^/?/2 + 03^%) 



02 



{o^iVm ~ oi3^) 



03 



/V3 

010203 

s/mmv3 



(aiy/rn - a2y/m) 



("ia/^ - a2^/m - "3 V^) • 



Remember that we assumed that ai^/rji > 
'^2v^i 013 V^' 'which implies that the relation ()58ap is 
violated. Thus, we have a > and 6 > 0. Therefore, the 
remaining two eigenvalues are also non-negative. Thus, 
xf are the solution of the relaxed two-state problem. 

If xf are all non-negative, xf is the optimal solution 
Xi of the original three-state problem. 



1 



1-1 



0:2 = 1 - 



^3 = 1 - 



(7^1(01102)1 + ^/^1(03101)1), 



(67a) 



(V^l(0l|02)|-V%l(02l03)l), 

(67b) 
(^^^l(03|0l)|-V^|(02|03)|), (67c) 



and the maximum success probability Pmax is given by 



Pn 



1-20^^1(01 I 02)1 +20^1(02 I 03)1 



2V^7i¥l(03|0i)|. 



(68) 




Note that apparent asymmetry of the signs of terms in 
the above results is due to our assumption that aiy/rji > 

If some of x^ are negative, we can still have analytical 
results by considering the appropriate two-state problem 
in the same way as in other cases. 



C. Examples 



We consider the set of the following three-state vectors: 



02} = 



COS (p3 sin 03 

sin 1^3 sin 9^ 

cos 6*3 

In the first example, the equal occurrence probabilities 
are assumed, and ip2 and (p^ are fixed to be 7r/3 and 7r/4, 
respectively. Figure [2] displays the optimal solution Xi 
and the maximum success probability as a function of 63 . 
As 03 varies from to 7r/2, the type of solutions changes. 
When 03 is small, all the three states have finite identi- 
fication probabilities Xi. As ^3 increases, the number of 
the states with nonvanishing Xi decreases. This rcfiects 
the fact that the state | (pz ) approaches the state | <t)2 ) 
and it becomes difficult to discriminate between these 
two states. When 9^ is in the range from 0.18(rad) to 
l.Ol(rad), xi and X2 are constant. That is when the 
three lengths ai^Jfjl satisfy the triangle condition ([58]) in 
the case of F > 0. In such a case, xi is expressed only by 
the mutual inner products as in Ea.(|55)). 



1 r- 



Xi 



X2 



cos <P2 COS (/33 
C0S(V32 - VsY 

cosy72Cos(y2 - Vz) 

cos (^3 



cos (p3 cos ((^2 -Vi) ■2a 

xz — ^ sm 9z- 

cos (^2 

It is clear that x\ and X2 are independent of ^3. Thus, 
x\ and X2 are constant in this range. 

In the second example, the occurrence probabilities 
are varied while the three states are fixed by setting 
ip2 = 7''/3, ip3 = 7r/4, and ^3 = 7r/5. Assuming 771 =772, 
we vary 773 from to 1 . Figure [3] displays the optimal 
solution Xi and the maximum success probability as a 
function of 773. For a sufficiently small 773, the optimal 
solution gives 2:3 = 0. As 773 approaches 1, only the state 
I ^3 ) survives as expected. However, the maximum suc- 
cess probability does not approach 1. This is because we 
assumed 77^ 7^ and imposed the no-error condition ^ 
for all states even if some of 774 are very small. One in- 
teresting thing is that if F > and the triangle condition 
(I55|) is satisfied, Eq. (|63p shows that Xi docs not depend 
on r/i! The optimal measurement is independent of the 
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FIG. 2: (Color online) Example 1. Xi (upper part) and the 
maximum success probability (lower part) vs. ^3. The angles 
</52 and ip3 are fixed to be 7r/3 and 7r/4, respectively. 



occurrence probabilities. This is the reason why all the 
Xi are constant when 773 ranges from 0.05 to 0.7. 

Before concluding the section, we mention that the 
problem of unambiguous discrimination can be numer- 
ically solved since the problem is formulated as semidef- 
inite programming. We performed a numerical analysis 
for the above examples using the semidefinite program- 
ming package SDPA [231, ^^id ^^ compared the two re- 
sults just to be sure. We verified that the numerical 
results reproduced our analytical ones within numerical 
errors. 



IV. DISCUSSION 

In the preceding section, we have obtained complete 
analytical solutions for unambiguous discrimination of 
three pure states with real mutual inner products and 
general occurrence probabilities. Here, we examine what 
difficulties arise when the inner products are complex. 
Though our method is not applicable to general complex 
cases, we will show that there is one instance for which 
an analytical result can be obtained. 

As in the preceding section, we define F in the following 
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FIG. 3: (Color online) Example 2. Si (upper part) and the 
maximum success probability (lower part) vs. 773. i]i — 772 
is assumed. The three states are fixed by setting (^2 ~ tt/S, 
fs = 7i"/4, and ^3 = tt/5. 



way: 



; 01 I 02 ) ( (/'2 I (/-a ) ( 03 I (/-i 



irle' 



3ie 



(69) 



where < < 27r/3. We can evidently choose phases 
of the three states so that the phases of inner products 
( 01 I 02 ) , ( 02 I 03 ) , and ( 03 I 01 ) are equal to the same 
value 9. Then, we can express the three inner products 
in the following way: 



1 102) = aia2e''', 


(70) 


2 103) = a2a3e''', 


(71) 


3 |0i) = aaaie''^, 


(72) 



where real positive numbers a^ arc now defined in terms 



of the modulus of the inner products, 



"1 



a2 



as 



n\<P2, 



^3\9l. 



1(02103)1 



' |(0l|02)||(02|03)| 
l(03|0l)l 



»2 I <?>3 , 



»3 h^l 



l(0l|02)l 



(73) 
(74) 
(75) 



Writing Yij = (yi\yj) as before, we can express d{^ 
tryiV) as follows: 

d=l + {{A\D\A)), (76) 

where | ^)) is the three-component vector defined as 



ail 2/1) 

1^))= I «2|2/2) 

aslvs) 

and D is the 3x3 matrix given by 

e-*" e' 



(77) 



D 



e- 



(78) 







The Hermitian matrix D is readily diagonalized as 

C/ti?C/ = diag(Ai,A2,A3), (79) 

by the unitary matrix given by 



U 




(80) 



where uj = e' 3 '^ is the cubic root of unity. The eigenvalues 
of D are found to be 

Ai = 2cos6', A2 = 2 cos (6* - 27r/3) , 

A3 = 2 cos (6* + 27r/3) . (81) 

We can now express d of Eq. ([75)) in the following form: 

(82) 
where the three vectors | at ) are defined by 



rf=l+g(Ai||ai)|VA2||a2)|%A3||a3)f 



\Z) I =V3U^\A)) 
las) 



ailyi) +0:21^2) +"312/3) 
= I "iIj/i) +i^^a2|2/2) +wa3|y3) 
aij j/i) +uja2\y2) +uj'^a3\y3) 



(83) 



Equation ([82)) shows that d is a linear combination of 
squared norms of the three vectors | a^ ) . Remember that 
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yi ) is arbitrary with only one constraint on its norm, 
\yi)\ = \/Wi- The maximum and minimum of the norm 
of each | ai ) can easily be obtained: The minimum of 
|| Oj;)! is if the three lengths ai^/rji satisfy the trian- 
gle condition (|58p . Otherwise, the minimum is given by 
ctiy/rfi-' a2y/rJ2 — ctSy/vE, where we assume aiy/rji is the 
greatest of all ai^Jrfi. The maximum of | a^ ) is evidently 
given by aiy/ryl + a2y^ + a^y/rj^,. We also notice that 
the sum of three squared norms of | a^ ) is constant. 



3 

El 



= 3((^|A))=3^ 



aJVi 



(84) 



Using this identity, we can rewrite d in terms of any pair 
of squared norms of | a^ ) , which is still difficult to ana- 
lytically minimize. However, if two of the eigenvalues Ai 
are equal, d can be expressed by the squared norm of a 
single vector | a^ ) , which can be minimized very easily. 
We find that this occurs if and only if = or = 7r/3. 
Clearly, these two cases correspond to the case of real F, 
which was solved in the preceding section. It can be ver- 
ified that we obtain the same results as before in these 
cases. 

There is one instance for which we can easily obtain an 
analytical result even if F is complex. This is the case in 
which three aiy/rjl are equal to each other. Eliminating 
1 1 ^3 ) I by using Eq. ([84|) . we have 



where D* is the complex conjugate of D, and the inequal- 
ity follows because A3 is also the smallest eigenvalue of 
D* . Therefore, the set of xf- given by Eq. (|88l) is the so- 
lution of the relaxed problem. If all x^ are non-negative, 
they are also the solution of the nonrelaxed problem. The 
maximum success probability is then given by dmin. 



l + 2cos(' 



2 

— TT, / 

3 'U 



3 

'^ aim 



(90) 



Otherwise, the problem is reduced to a certain two-state 
problem. 

This result is a generalization of the symmetric states 
considered in Ref. |ll| . though the number of states is 
limited to three here. For the symmetric states, the ab- 
solute values of all mutual inner products are the same, 
implying a constant a^, and the occurrence probabilities 
r]i are assumed equal. In this case, xi' are shown to 
be the same and nonnegative, and the maximal success 
probability is given by 



Pmax = 1 + 27 C0S(6' + -Tt), (91) 

where 7 = |(0i | 02 )| = |( (/>2 | '?!>3 )| = |( 03 I 0i )|- 



V. SUMMARY AND CONCLUDING REMARKS 



d = 1 + A3 ^ aim 

1=1 

+ ^((Ai-A3)||ai)|' + (A2-A3)||a2)|'). (85) 

We find that the two norms | ai ) and | 02 ) can si- 
multaneously be if aiy/fji is constant. This is possible 
by choosing \yi), for example, to be single-component 
complex vectors so that three a^ | y^ ) make an equilateral 
triangle in the complex plane. 



1 2/1 > = Vm, 12/2) = ^/?/2w^ 12/3) = Vm^- 



(86) 



Note that Ai > A3 and A2 > A3. Thus, the minimum of 
d is given by 



dmin = 1 + 2 C0S( 



;-)E 



aim- 



(87) 



i=l 



The attainability condition (|26| requires that 

2 

xf = l + 2a1cos{e + -?:). 
o 

The positivity of A^— X^ can be verified in the following 
way: 

N- X^ 

= diag(ai, ^2, as) (D* - A3I) diag(ai, a2, as) 
> 0, (89) 



We have established some general properties of un- 
ambiguous discrimination of n linearly independent pure 
states. By formulating the problem as semidcfinitc pro- 
gramming, we have shown that the solution of the prob- 
lem is unique. In unambiguous discrimination, the opti- 
mal strategy may produce a vanishing probability of iden- 
tifying the input state with some of the states. To deal 
with this possibility, it has been shown to be very con- 
venient to consider a problem in which some constraints 
on the variables are relaxed. Applying our method to 
the case of three pure states, we have obtained complete 
analytic solutions if their mutual inner products are real. 
Unfortunately, our method is not generally applicable to 
the case of complex inner products. However, we have 
found a class of ensemble of three pure states with com- 
plex inner products for which we obtain an analytical 
solution. 

It will also be of interest in future studies to extend 
our method to discrimination of three pure states with 
error margin [9|,|l0| and to discrimination of three unitary 
processes |3g |. 



Appendix: Uniqueness of solutions in the n-state 
problem 

We will show the solution of the n-state problem is 
uniquely determined. The following two propositions 
hold: 
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Proposition 1. If x and x' are both solutions of the 
n-state problem, then Xig = x'- for at least one Iq. 

Proposition 2. Consider the n-state problem with an 
additional condition: One of the components Xi^ is fixed 
to be some value. We assume a solution exists. Then, the 
remaining variables are either or uniquely determined 
by a certain r -state problem, where 1 < r < n — 1. 

With these two propositions, the uniqueness of the so- 
lution of tlie n-state problem follows by induetion on the 
number n of the states. It is evident when n = 1. We 
assume that the solution is unique in the r-state problem 
for r < n — 1. Let x and x' be two solutions of the n- 
state problem. By proposition (I), we have Xi^ — x[ for 
a certain jq- Proposition (II) states that the remaining 
variables are determined by a certain r-state problem, 
where 1 < r < n — 1. By the assumption, we can con- 
clude X = x' . 

Now we present the proof of the two propositions. 

Proof of Proposition 1. li N — X > 0, then we can add 
a sufficiently small positive quantity to any component 
Xi so the condition A'' — X > is still satisfied while the 
success probability increases. Therefore, det{N — X) = 
if a; is a solution. Since a convex linear combination of 
the two solutions is also a solution, we have 

dct[N - kX - {1 - k)X'] = (0<K<1). (A.l) 

The coefficient of k" of this determinant is given by 

{x[~Xi){x2-X2)---{x'„-Xri), (A.2) 

which should be 0, implying Xi^ ~ x[ for a certain zq. I 

Proof of Proposition 2. Without loss of generality, we as- 
sume zo = n. We write the condition d?]) as 



n-l 

A-^a;,|0O('^.|>O, 

i = l 
A = ly - Xn\(i>n){4'n \- 



(A.3) 

(A.4) 



The assumption that a solution exists to the problem 
implies A > 0. We first consider the case A > 0. The 
condition (|A.3|) can be written as 



lv-^x,A-i/2|0,)(0»|A"'/'>O. 



(A.5) 



The states {A~^/^| (jii )}^Ii are linearly independent. By 
V , we denote the in — l)-dimensional subspace spanned 
by them. We define the following two sets of n — 1 states 
in^': 



\4>[) = 



Q\' 



v/(^7|Q^ 



(i = l,...,n- 1), 



l0n = v/('/'dQI0.>A-i/'|^.) (,^i^...^„_ 



(A.6) 

1), 
(A.7) 



where Q is the projector onto the subspace V . It can 
be readily verified that {| (/>^ llfJi^ and {| ^[ )}^r/ are re- 
ciprocal in the sense that ( 0^ | 0'- ) — Sij . The condition 
(IA.5I) takes the form given by 



1V,-J2^^\4>^){4>'^\>0, 



(A., 



where x'j = Xi/{4>i\Q\4>i). The success probability to be 
maximized can be expressed as 

n n — 1 

p^^riiXi^rinXn + N^ri-x'i, (A. 9) 



i=l 



where 



ijl = n,{cj,,\Q\^,)/N, 



Tl-1 



N^J2^^ 



Thus, maximization of p in Eq. (|A.9p subject to the con- 
ditions (jA.Sp and x^ > is equivalent to the (n— l)-state 
problem for the ensemble 1??^, | 0^ )}"Ji^- 

Next, we consider the case in which A is singular. We 
observe that this happens if and only iixn ~ l/{(f>n\'Pn), 
and that A is the projector such that 

A|0,) = |0*> (i = l,...,n-l), (A.IO) 

A|<^„)=0. (A.ll) 

We also find that Xi should be unless {(l>n\4'i) = 0- 
Without loss of generality, we assume that 







nonzero i 



1,- 



l,...,n- 1, 



(A.12) 



where < r < n—1. If r = then all Xi are 0. We assume 
r > 1. Let T^' be the subspace spanned by the linearly 
independent states {| (p'^ )}l=i- Since 1^' is evidently con- 
tained in the support of A, the condition (|A.3[) can be 
written as 



^a;,|0,)(<^»l>O- 



(A.13) 



4=1 



The states which are reciprocal to the set {| (pi )}[^i are 
not generally normalized. However, it is always possible 
to obtain the set of normalized states {| 0i)}[^i which 
is reciprocal to the set {ci\ <j>i )}[^i by choosing the coef- 
ficients Ci appropriately. Note that the resultant states 
I 4)[ ) are unique up to a phase factor. Thus, maximizing 
the success probability subject to the conditions (|A.13|) 
and Xi > is equivalent to the r-state problem for the 
states {| (p'i )}l^i with properly modified occurrence prob- 
abilities ril. ■ 
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